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Toxicity  of  50  organic  chemicals  to  microorganisms  was  determined  using  the 
respirometeric  approach.  Using  this  experimental  database,  models  for  predicting 
toxicity  (IC50  values)  were  developed  using  QSAR  techniques.  Toxicity  measurements 
were  also  made  for  ten  binary  mixtures,  and  sixteen  multi-component  mixtures.  The 
joint  effects  of  organic  chemicals  in  mixtures  were  analyzed  by  three  different 
approaches.  Using  the  QSAR  model  developed  from  single  chemical  studies,  an  approach 
was  developed  to  analyze  and  predict  joint  effects  of  chemicals  in  mixtures.  The 
results  of  this  study  indicate  that  the  joint  effects  could  be  considered  simply 
additive  for  the  different  classes  of  chemicals  tested. 
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Abstract 

Toxicity  of  50  organic  chemicals  to  microorganisms  was  determined  using  the 
respirometeric  approach.  Using  this  experimental  database,  models  for  predicting 
toxicity  (IC5Q  values)  were  developed  using  QSAR  techniques.  Toxicity  measurements 
were  also  made  for  ten  binary  mixtures,  and  sixteen  multi-component  mixtures.  The 
joint  effects  of  organic  chemicals  in  mixtures  were  analyzed  by  three  different 
approaches.  Using  the  QSAR  model  developed  from  single  chemical  studies,  an 
approach  was  developed  to  analyze  and  predict  joint  effects  of  chemicals  in  mixtures. 
The  results  of  this  study  indicate  that  the  joint  effects  could  be  considered  simply 
additive  for  the  different  classes  of  chemicals  tested. 

Introduction 

Acute  and  chronic  toxicity  testing  is  a  major  component  in  the  NPDES  permitting 
process.  The  concept  of  whole  effluent  toxicity  testing  has  been  introduced  into  this 
program  due  to  the  realization  that  a  mixture  of  several  chemicals  may  exhibit  greater 
toxicity  than  they  would  individually.  While  current  Water  Quality  Standards  are  based 
on  single  chemical  toxicity  assays,  in  future,  controls  may  be  set  based  on  the  joint 
effects  of  mixtures  of  two  or  more  chemicals.  Non-point  sources,  industrial  effluents, 
leachates  and  contaminated  groundwaters  are  all  known  to  contain  several  chemicals 
in  mixtures.  Thus  an  ability  to  analyze  and  predict  joint  effects  of  mixtures  of  chemicals 
on  microorganisms  and  other  aquatic  life  forms  will  be  of  considerable  benefit  in 
managing  the  environmental  hazards  of  synthetic  chemicals. 

Several  ecological  researchers,  notably  from  Europe,  have  studied  the  effects  of 
mixtures  of  chemicals  on  fish  (Ref  1  -  20).  Hardly  any  studies  have  been  reported  on 
the  joint  effects  of  several  chemicals  on  microorganisms.  As  microorganisms  are 
employed  in  municipal  waste  treatment  by  environmental  engineers,  and  are  also 
present  in  the  natural  environment,  it  would  be  of  interest  to  be  able  predict  such 
effects  on  microorganisms.  A  study  is  underway  in  our  laboratories  to  determine  and 
predict  joint  effects  of  binary  and  multiple  chemical  mixtures  on  three  classes  of 
organisms  of  interest  to  environmental  engineers.  This  interim  report  covers  the  results 
of  the  Phase  I  study  during  the  first  year. 
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Objectives  of  Phase  I  Study 

The  ultimate  objective  of  this  3-year  research  is  to  develop  an  approach  to  predict 
the  joint  toxic  effects  of  mixtures  of  organic  chemicals  to  microorganisms.  Towards  this 
end,  the  following  tasks  were  identified  for  the  first  year: 

a)  measure  single  chemical  toxicity  to  a  surrogate  culture  of  microorganisms; 

b)  develop  a  QSAR  model  to  predict  single  chemical  toxicity  to  the  surrogate 
organisms; 

c)  measure  and  analyze  toxicity  of  binary  mixtures  to  surrogate  organisms  and 
verify  simple  additivity; 

d)  measure  and  analyze  toxicity  of  multi-component  mixtures  to  surrogate 
organisms  and  verify  simple  additivity; 

e)  develop  an  approach  to  predict  joint  effects  of  mixtures  of  organic  chemicals 
to  the  surrogate  culture  based  purely  on  the  molecular  structures  of  the 
components  of  the  mixtures. 

Work  during  the  second  and  third  years  will  focus  on  the  effects  of  mixtures  on 
activated  sludge  and  anaerobic  cultures  with  a  view  to  predict  joint  effects  on  them 
using  molecular  structures  of  the  components  and  the  surrogate  culture  results. 

Several  techniques  and  approaches  in  toxicity  modeling  reported  in  the  literature 
as  individual  studies  are  integrated  in  this  study,  ^me  of  the  features  of  this  study  are 
use  of  a  20-reactor,  computer  integrated  respirometer  for  toxicity  determination;  use  of 
the  probit  transformation  method  to  determine  IC50;  and  use  of  QSAR  approach  to 
predict  joint  effects  of  binary  mixtures  and  multi-component  mixtures. 

This  interim  report  covers  the  work  done  in  Phase  I  during  the  first  year.  Results 
reported  below  include; 

a)  developing  an  experimental  procedure  for  IC5Q  determination  using  the 
respirometeric  approach; 

b)  evaluation  of  the  use  of  Polytox  (poibac  Corp.)  cultures  as  surrogate 
microorganisms  in  rapid  toxicity  determination; 

c)  development  of  QSAR  model  for  single  chemical  toxicity; 

d)  testing  the  hypothesis  that  joint  effects  of  the  chosen  chemicals  are  simply 
additive  in  binary  and  multi-component  mixtures; 

e)  establishment  of  a  methodology  to  analyze  joint  effects  in  binary  mixtures 
and  multi-component  mixtures; 

f)  development  and  validation  of  a  QSAR-based  approach  to  predict  the 
concentrations  of  the  components  in  mixture  that  would  cause  50% 
inhibition. 
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Experimental  Approach 

All  tests  were  conducted  using  research  grade  chemicals  as  supplied  by  the 
manufacturers  without  any  further  purification.  The  toxicity  tests  were  run  on  a  20- 
reactor  computer  interfaced  N-Con  Respirometer  Model  as  detailed  in  Appendix  I.  The 
test  procedure  is  detailed  in  Appendix  II.  The  percent  inhibition  caused  by  a  toxicant  at 
a  given  concentration  was  determined  by  comparing  the  oxygen  uptake  rate  of  a 
toxicant-free  control  reactor  against  the  rates  of  eight  other  reactors  spiked  with 
different  concentrations  of  the  toxicant.  This  rate  was  in  turn  obtained  from  the  slopes 
of  the  linear  portion  of  the  oxygen  uptake  curves  generated  by  the  respirometer  for 
each  reactor.  The  %  inhibition  values  were  then  transformed  to  probit  values 
according  to  Finney  (1964).  From  these  probit  vs.  concentration  plots,  the  ICsq  values 
were  then  determined  [IC50  is  the  concentration  of  the  toxicant  at  which  a  50% 
inhibition  is  caused].  The  above  procedure  is  illustrated  schematically  in  Fig  1 . 
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Figure  1 .  Schematic  Illustration  of  Procedure  to  Find  IC5P 


In  the  single  chemical  studies,  50  chemicals  belonging  to  several  different 
congeneric  classes  with  a  diversity  of  molecular  structures  were  chosen  for  testing. 
These  chemicals  spanned  a  wide  range  of  aqueous  solubility,  Henry’s  Constant  and 
octanol-water  partition  coefficient.  Such  a  variety  of  chemicals  would  enhance  the 
robustness  and  the  utility  of  the  QSAR  models. 

In  the  binary  mixture  tests,  n-octanol  was  chosen  as  one  of  the  toxicants  and  the 
following  chemicals  as  the  second  component:  chlorobenzene,  2,4  dimethyl  phenol, 
1,1,1  trichloroethane,  cyclohexane,  n-butyl  amine,  n-amyl  acetate,  ethanol, 
chlorodibromomethane,  cyclohexanone  and  ethanolamine.  For  each  combination  of 
the  two  components,  a  matrix  of  four  octanol  concentrations  by  four  second 
contaminant  concentrations  was  assayed.  The  four  test  concentrations  were  chosen 
as  0.4, 0.6, 0.8  and  0.9  of  the  individual  IC50  concentrations. 
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In  the  multi-component  mixture  tests,  ten  mixtures  each  containing  10  different 
chemicals  mixed  in  equitoxic  proportions  were  assayed.  For  each  mixture,  one  control 
reactor  and  six  spiked  reactors  were  run.  The  six  reactors  received  0.04,  0.06,  0.08, 
0.1,  0.12  and  0.14  Toxic  Units  of  each  of  the  ten  components.  Additional  six  testing 
mixtures,  each  containing  8  different  chemicals,  were  also  assayed.  In  this  case,  the 
six  spiked  reactors  received  0.05,  0.075,  0.1,  0.125,  0.150  and  0.175  Toxic  Units  of 
each  of  the  components. 

Modeling  Approach 

The  single  chemical  toxicity  results  from  the  50  chemicals  were  used  to  develop  a 
QSAR  model.  Molecular  connectivity  indexes  were  calculated  for  the  50  chemicals 
following  the  algorithms  developed  by  Kier  and  Hall  and  modified  by  Nirmalakhandan 
(1988).  Simple  and  multiple  step-wise  regression  analysis  procedures  were  used  to 
derive  the  QSAR  model  with  IC5Q  values  as  the  dependent  variable.  The  IC5Q  values 
calculated  from  the  QSAR  models  were  then  compared  with  the  experimentally 
measured  values. 

In  the  binary  mixture  studies,  isobolograms  were  first  constructed  to  test  the 
hypothesis  that  the  joint  effects  were  according  to  the  simple  additivity  model.  Then, 
the  QSAR  models  developed  from  single  chemical  was  used  to  predict  the  joint  toxicity 
caused  by  the  second  chemical  in  the  presence  of  n-octanol.  The  predicted  values 
were  compared  with  the  experimental  results. 

In  the  multi-component  mixture  studies,  the  joint  effects  were  analyzed  using  three 
concepts:  Toxic  Unit.  Additivity  Index,  and  Mixture  Toxicity  Index.  The  validity  of  these 
concepts  was  further  verified  using  the  results  of  the  8-component  testing  set.  Finally, 
the  QSAR  models  developed  from  single  chemical  tests  were  used  to  predict  the 
concentrations  of  the  components  in  the  8-component  mixture  that  would  cause  50% 
inhibition.  These  predicted  concentrations  were  then  compared  with  the 
experimentally  measured  concentrations. 


5 


Results  and  Discussion 
a)  Single  chemical  results* 

The  50  chemicals  used  in  this  study  are  listed  in  Table  I.  This  Table  also  contains 
the  respective  IC50  values  and  the  corresponding  r^  found  from  the  probit  vs. 
concentration  plots  from  which  the  IC50  values  were  determined.  These  r^  values 

range  from  0.7  to  0.98,  indicating  a  reasonably  good  database.  The  distribution  of  r^  is 
shown  in  Fig  2.  On  preliminary  evaluation,  the  following  four  chemicals  were  found  to 
be  "outliers”:  1,4  dichlorobenzene.  1.2  dichloroethylene.  triethanolamine  and 
bromodichloromethane.  At  this  point,  no  explanation  could  be  postulated  for  the 
outlying  tendency  of  these  chemicals.  Other  researchers  have  found  similar 
peculiarities  too.  These  chemicals  were  therefore  excluded  from  the  analysis  and  the 
remaining  46  were  used  in  the  OSAR  model  development. 

The  first  order  molecular  connectivity  index,  was  found  to  fit  the  individual  IC5Q 
values  (in  mM/L)  for  the  alkanes,  aromatic  compounds,  alcohols,  esters,  ketones, 
amines: 

Alcohols,  ketones  and  esters: 

log  IC50  »  3.690  -  0.896  ^x'' 

n  =  14;  r  =  0.954;  =  0.910;  RMS  residual  =  0.246. 

Alkanes: 

log  IC5o  =  1.851  -  0.765  ^X'' 

n  =  5;  r  =  0.999;  r^  =  0.999;  RMS  residual  =  0.018. 

Amines  and  acids: 

log  IC5o  =  1.045 -0.470  ■'x'' 

n  =  6;  r  =  0.957;  r^  =  0.915;  RMS  residual  «  0.101. 

Aromatics: 

log  IC50*  3.258- 1.133  ^X'' 
n  =  9;  r  *  0.852;  r^  s  0.726;  RMS  residual  *  0.31 1 . 

In  the  case  of  halogenated  aliphatics,  the  zero  order  molecular  connectivity  index, 
^X'^.  fitted  the  data  set  more  closely: 

Halogenated  aliphatics: 

log  IC50  =  2.670  -  0.448  ^x'' 
n  «  12;  r  =  0.942;  ■  0.887;  RMS  residual  »  0.141 . 
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Table  I.  Chemicals  Assayed  and  Single  Chemical  Results 


ID  Chemical 

# _ 

1.  Benzene 

2.  Toluene 

3.  Xylene 

4.  Ethylbenzene 

5.  Chlorobenzene 

6.  1,2  Dichlorobenzene 

7.  1,3  Dichlorobenzene 

8.  1,4  Dichlorobenzene 

9.  1,2,4  Trichlorobenzene 

1 0.  2,4  Dimethyl  phenol 

11.  Methylene  chloride 

12.  Dibromomethane 

13.  Carbon  tetrachloride 

14.  1,2  Dichloroethane 

15.  1,1,1  Trichloroethane 

16.  1,1,2,2  Tetrachloroethane 

17.  1,2  Dichloropropane 

18.  Bromochloromethane 

1 9.  Bromodichloromethane 

20.  Chlorodibromomethane 

21.  Ethylene  dibromide 
22. 1,2  Dichloroethylene 

23.  Trichloroethylene 

24.  Tetrachloroethylene 

25.  Cyclohexane 

26.  Pentane 

27.  Hexane 

28.  Heptane 

29.  Octane 

30.  Bis  (2*chloroethyl)  ether 

31.  Ethanol 

32.  Propanol 

33.  Pentanol 

34.  Octanol 

35.  n-Butyl  acetate _ 


Type 

IC50 

[mg/Ll 

m 

ARO 

685 

0.898 

ARO 

207 

0.870 

ARO 

140 

0.851 

ARO 

220 

0.950 

ARO 

350 

0.966 

ARO 

135 

0.939 

ARO 

40 

0.903 

ARO 

6 

0.908 

ARO 

23 

0.810 

ARO 

240 

0.953 

HAL 

1,750 

0.845 

HAL 

1,110 

0.848 

HAL 

325 

0.701 

HAL 

685 

0.979 

HAL 

415 

0.850 

HAL 

180 

0.949 

HAL 

500 

0.954 

HAL 

1,800 

0.948 

HAL 

90 

0.916 

HAL 

425 

0.959 

HAL 

520 

0.939 

HAL 

350 

0.783 

HAL 

500 

0.979 

HAL 

175 

0.923 

ALK 

74 

0.964 

ALK 

70 

0.793 

ALK 

38 

0.751 

ALK 

18 

0.739 

ALK 

8 

0.857 

AKE 

1,600 

0.886 

AKE 

40,000 

0.964 

AKE 

7,200 

0.942 

AKE 

2,325 

0.821 

AKE 

126 

0.968 

AKE 

3,750 

0.923 

Log  IC50  [mM/U 
Exp.  Calc.  Error 


0.94 

0.99 

-0.05 

0.35 

0.53 

-0.18 

0.12 

0.06 

0.06 

0.32 

-0.11 

0.43 

0.50 

0.46 

0.04 

-0.04 

-0.10 

0.06 

-0.57 

-0.09 

-0.48 

-1.39 

7 

-0.90 

-0.63 

-0.27  * 

0.29 

-0.10 

0.39 

1.31 

1.34 

-0.03 

0.81 

0.60 

0.21  * 

0.32 

0.42 

-0.09 

0.84 

1.02 

-0.18 

0.49 

0.48 

0.02 

0.03 

0.13 

-0.10 

0.65 

0.63 

0.01 

1.15 

0.97 

0.18 

-0.26 

7 

0.31 

0.15 

0.17 

0.45 

0.37 

0.08 

0.56 

7 

0.58 

0.67 

-0.09 

0.02 

0.19 

-0.17 

-0.06 

-0.06 

0.01 

-0.01 

0.01 

-0.02 

-0.36 

-0.38 

0.02  • 

-0.75 

-0.76 

0.01 

-1.15 

-1.14 

-0.01 

1.05 

0.84 

0.21  * 

2.94 

2.79 

0.15 

2.08 

2.33 

-0.25  * 

1.42 

1.41 

0.01 

-0.01 

0.04 

-0.05 

1.51 

1.07 

0.44  ** 

Code  for  type:  ARO>  aromatics;  HAL-  halogenated  aliphatics;  ALK-  alkanes: 
AKE-  alcohols,ketones  and  esters;  AMI-  amines. 

Note:  *  error  >1  std.  dev.;  **  error  >  2  std.  dev.;  Z  excluded  from  analysis. 


Table  I  (contd.) 


ID  Chemical 
# 

Type 

IC50 

[mg/L] 

rA2  __ 

Log  IC50  fmM/L] 

Exp.  Calc.  Error 

36.  Isobutyl  acetate 

AKE 

1,600 

0.955 

1.14 

1.20 

-0.06 

37.  n-Amyl  acetate 

AKE 

440 

0.980 

0.53 

0.61 

-0.08 

38.  Ethyl  acetate 

AKE 

5,400 

0.916 

1.79 

1.98 

-0.19 

39.  Acetone 

AKE 

48,000 

0.864 

2.92 

2.62 

0.30  * 

40.  Methyl  ethyl  ketone 

AKE 

1,900 

0.885 

1.41 

1.90 

-0.49  ** 

41.  Methyl  isobutyl  ketone 

AKE 

2,600 

0.905 

1.42 

1.32 

0.09 

42.  Methyl  n-propyl  ketone 

AKE 

4,500 

0.883 

1.72 

1.65 

0.07 

43.  Cyclohexanone 

AKE 

3,750 

0.964 

1.58 

1.51 

0.07 

44.  n*Butyl  amine 

AMI 

90 

0.952 

0.09 

0.05 

0.04 

45.  t-Butyl  amine 

AMI 

85 

0.919 

0.07 

0.21 

-0.14 

46.  Diethylamine 

AMI 

104 

0.769 

0.15 

0.05 

0.11 

47.  Acetic  acid 

AMI 

287 

0.744 

0.68 

0.61 

0.07 

48.  Cyclohexylamine 

AMI 

60 

0.885 

*0.22 

-0.20 

-0.02 

49.  Ethanolamine 

AMI 

160 

0.955 

0.42 

0.47 

-0.05 

50.  Triethanolamine 

AMI 

900 

0.829 

0.78 

? 

Average  0.893 


Code  for  type:  ARO-  aromatics;  HAL*  halogenated  aliphatics;  ALK*  alkanes; 
AKE*  alcohols, ketones  and  esters;  AMI*  amines. 

Note:  *  error  >1  std.  dev.;  **  error  >  2  std.  dev.;  ?  excluded  from  analysis. 


Coefficient  of  Determination,  r  ^2  Predictive  Error 


Figure  2.  Frequency  Distribution  of 
Coefficient  Determination  and  Predictive  Error  in  Table  I 
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Detailed  statistical  parameters  of  these  equations  are  included  in  Appendix  III. 
These  models  explain  about  90%  of  the  variance  in  the  experimental  data  for  alcohols, 
ketones,  esters,  alkanes,  amines  and  aromatics  and  over  70%  in  the  case  of  the 
halogenated  aliphatics.  The  average  RMS  residuals  of  0.163  suggests  that  these 
equations  can  predict  IC5Q  within  a  factor  of  1 .5.  Except  in  the  case  of  the  aromatic 
compounds,  this  degree  of  precision  is  comparable  to  the  uncertainty  of  the 
experimental  IC5Q  values.  The  IC5Q  values  calculated  by  these  equations  are 
tabulated  in  Table  I  and  compared  against  the  respective  experimental  values  in  Fig  3, 
from  which  the  overall  quality  of  the  agreement  between  the  two  can  be  seen  to  be 
comparable  to  the  experimental  uncertainty  in  the  IC50  values,  with  r^  =  0.954.  The 
distribution  of  the  predictive  error  is  shown  in  Fig  2. 


Calculated  log  IC50 


Figure  3.  Comparison  Between  Experimental  and  Calculated  log  IC5Q  Values  [mM/L]. 
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b)  Binary  mixture  results- 

The  isobologram  developed  from  the  experimental  IC5Q  of  the  binary  mixtures  is 
shown  in  Fig  4.  The  experimental  data  points  for  the  six  toxicants  falling  on  or  around 
the  diagonal  confirms  the  hypothesis  of  simple  additivity.  Additional  tests  are  being 
done  now  to  check  these  results  and  include  more  chemicals. 


Toxic  Units  of  Octanol 


Figure  4.  Isobologram-  Octanol  vs.  Six  Other  Chemicals 
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c)  Multi-component  mixture  results- 

The  results  of  the  multi-component  tests  are  presented  in  graphical  form  in 
Appendix  IV-a  for  the  10-component  mixtures  and  in  Appendix  IV  b  for  the  8- 
component  mixtures.  The  r^  in  these  plots  range  from  0.860  to  0.995,  indicating  a 
satisfactory  database.  The  IC5Q  values  obtained  from  these  plots  were  used  to  analyze 
the  joint  effects  by  the  Toxic  Units,  Additivity  Index  and  Mixture  Toxicity  Index  concepts 
as  summarized  in  Table  III. 


Table  III.  Analysis  of  Joint  Effects  of  Multi-component  Mixtures 


Mixture 

N» 

Chemicals  in  Mixture 

Joint  Effects  Analyzed  bv  1 

ID  #  of  components 

Total 

n 

Additivity  Index 
AI.M-1 

Toxicity  Index 
MTI  =  1  -  loqM/loon 

IOC-1 

4.5,10.36,32,33,12,18.1.2 

10 

0.10 

0.96 

IOC-2 

4,5.10.36.32.33.12,18.22.23 

10 

1.00 

0.00 

1.00 

IOC-3 

40.41.35.36.32.33.4.5.10.17 

10 

1.14 

0.14 

0.94 

IOC-4 

40.41,35.36,32.33,4,5.10.2 

10 

0.98 

IOC-5 

40.41.35,36.32,33,31.43,34.17 

10 

1.31 

0.31 

0.88 

IOC-6 

40,41,35,36.31.43.12.18.1.2 

10 

1.10 

0.96 

IOC-7 

40,41,43.31.32,33.12.18,22.23 

10 

1.15 

0.15 

0.94 

IOC-8 

40.41,35,36,4.5.17.43,34,17 

10 

0.90 

•0.10 

1.05 

IOC-9 

40,41.35.36,4.5.17.18.1.2 

10 

0.90 

-0.10 

1.05 

10C-10 

40,41,43,4,5.17,12.18,22,23 

10 

0.90 

-0.10 

1.05 

8C-1 

40,41,35,36,32.33,12,30 

8 

1.19 

0.19 

0.92 

8C-2 

40,41,35,36,12,18,1,2 

8 

1.02 

0.02 

0.99 

8C-3 

40,41,35,36,32.33,22,23 

8 

1.07 

0.07 

0.97 

8C-4 

4,36.32.33,12.18.34,17 

8 

1.02 

0.02 

0.99 

8C-5 

4,10,36,32,33,18,22,23 

8 

0.92 

-0.08 

1.04 

8C-6 

40.35.21.15.4.5.10.2 

8 

0.93 

-0.07 

1.03 

Averaae 

1.04 

0.04 

0.981 

Note:  ID  #  refer  to  chemicals  listed  in  Table  I; 

Simple  additivity  implied  by  M  s  1 ,  Al  =  0  and  MTI  =  1 . 

The  results  of  these  analyses  presented  in  Table  III  support  the  hypothesis  of 
simple  additivity  for  the  chemicals  assayed,  with  average  values  of  £TU  =  1.04±0.12, 
Al  =  0.04±0.12,  and  MTI  =  0.98±0.05. 
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The  results  from  the  8-component  multi-mixture  studies  were  used  to  test  the 
hypothesis  that  the  concentrations  of  the  components  of  a  mixture  causing  50% 
inhibition  could  be  predicted  from  single  chemical  results.  Two  different  approaches 
were  used  to  predict  the  concentrations  of  the  components  at  which  50%  inhibition 
would  be  caused: 

-  in  one  approach,  perfect  simple  additivity  (i.e  ZTU  =  1.0)  was  assumed  and  the 

corresponding  concentrations  of  the  components  in  the  mixtures  were 
determined  from  experimentally  measured  single  chemical  IC^q  values. 

-  in  the  second  approach,  perfect  simple  additivity  (i.e  ITU  =  1.0)was  assumed  and 

the  corresponding  concentrations  of  the  components  in  the  mixtures  were 
determined  from  QSAR  models  for  single  chemical  IC^q  values. 

Results  of  these  calculations  and  comparisons  are  presented  in  detail  in  Appendix 
V.  The  concentrations  predicted  by  these  two  approaches  are  compared  against  the 
actual  experimental  concentrations  in  Fig  5  for  the  48  sets  of  data  points  for  each 
approach  (six  different  mixtures,  each  containing  8  components).  The  concentrations 
calculated  by  the  first  approach  is  in  near  perfect  agreement  (r^  =  0.995)  with  the 
experimentally  measured  concentrations.  This  confirms  that  the  joint  action  is 
according  to  simple  additivity  mechanism  for  the  components  tested. 

In  the  second  approach,  the  agreement  appears  to  bo  “poor"  (r^*  0.799).  This 
deviation  stems  mainly  from  the  inadequacy  of  the  QSAR  models  and  partly  from  the 
slight  deviation  from  simple  additivity.  In  addition,  two  of  these  mixtures.  Mixture  3 
and  5,  contain  1 ,2  dichloroethylene  which  was  found  to  be  an  “outlier”  in  the  single 
chemical  studies  and  excluded  in  the  QSAR  model  development.  Thus  the  single 
chemical  IC5Q  predicted  by  the  QSAR  model  for  this  chemical  results  in  a  higher  error. 

Nevertheless,  it  must  be  appreciated  that  the  second  approach  does  not  require 
any  experimental  inputs  whatsoever  and  the  predictive  errors  are  within  a  factor  of  3, 
and  the  predictions  are  practically  of  the  same  order  of  magnitude  as  the  experimental 
values,  in  the  span  of  2  log  units  of  IC50  concentrations.  This  degree  of  precision  is 
comparable  to  the  work  reported  in  the  literature  for  single  chemical  studies,  and  can 
be  considered  sufficiently  adequate  for  predicting  toxicity  of  multi-component  mixtures. 
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0  12  3  4 

Experimental  Cone,  of  Components 
in  Mixture  Causing  50%  Inhibition  [log(mg/L)] 


Figure  5.  Comparison  Between  Experimental  and  Predicted  Concentrations 


Conclusions 

The  respirometeric  technique  for  determining  toxicity  data  has  been  demonstrated 
to  be  an  efficient  and  reliable  one.  Joint  effects  of  the  range  of  chemicals  assayed  can 
be  adequately  quantified  by  the  simple  additivty  concept  using  toxic  units.  The  results 
of  the  above  studies  show  clearly  that  OSAR  approaches  can  be  used  to  predict  joint 
toxic  effects  of  multiple  chemicals  on  microorganisms.  The  quality  of  the  predictions  by 
the  QSAR  models  for  the  range  of  chemicals  tested  in  this  study  are  comparable  to  the 
experimental  uncertainty  in  the  IC5Q  data. 

Additional  chemicals  belonging  to  a  wider  range  of  molecular  structures  need  to  be 
assayed  and  used  for  checking  the  validity  and  utility  value  of  these  models.  In 
addition,  similar  studies  on  activated  sludge  and  anaerobic  cultures  need  to  be  done 
to  establish  correlation  between  the  Polytox  surrogate  cultures  so  that  far  reaching 
predictions  can  be  made  to  maximize  the  benefits  of  this  research. 
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APDflndiy  l- 

Details  of  Respirometer  System 


The  respirometer  being  used  in  our  research  was  developed  here  at  New  Mexico 
State  University,  and  is  commercially  marketed  by  N-CON  Corporation,  Inc.,  NY.  The 
system  has  been  recently  modified  in  our  laboratory  to  work  either  in  the  aerobic  or 
anaerobic  mode.  The  reactors  in  this  system  are  maintained  at  constant  temperature 
and  pressure.  Changes  in  headspace  pressure,  due  to  gas  production  (or 
consumption),  are  sensed  by  a  pressure  (or  vacuum)  switch  and  are  converted  to  gas 
volume  using  Ideal  Gas  Laws,  reactor  volume,  temperature  and  type  of  gas  being 
exchanged,  and,  monitored  on  a  real  time  basis.  These  volumes  can  then  be  easily 
related  to  biological  activity  in  the  reactor.  A  brief  description  of  the  system  in  the 
aerobic  mode  is  as  follows. 

In  this  mode,  the  CO2  produced  is  absorbed  by  KOH  pellets  placed  in  the 
headspace.  Thus,  consumption  of  O2  results  in  a  vacuum  in  the  headspace.  A  vacuum 
switch  has  its  vacuum  side  connected  to  the  headspace.  The  pressure  side  of  the 
switch  is  connected  to  a  closed,  constant  pressure  tank,  thus  providing  a  steady 
reference  pressure,  eliminating  any  fluctuations  due  to  barometric/atmospheric 
variations.  When  the  pressure  differential  across  the  switch  exceeds  2.5  mm  H2O,  a 
signal  is  sent  through  the  data  acquisition  system  to  the  computer  and,  a  precise  pulse 
of  oxygen  from  an  oxygen  cylinder  is  injected  into  the  headspace.  The  computer  keeps 
track  of  the  number  of  pulses  (or  the  amount)  of  oxygen  supply  as  a  function  of  time. 
From  this  data,  oxygen  utilization  rate  can  then  be  established.  A  schematic 
arrangement  of  this  system  is  shown  below: 


Aopgndijs  H; 


In  this  study,  a  commercial  culture-  Polytox,  was  used  to  determine  IC50  values  of 
50  organic  chemicals.  One  vial  of  8  gr  of  Polytox  freeze  dried  cultures  was  dispersed 
in  250  mL  of  buffered  dilution  water,  prepared  according  to  Standard  Methods.  This 
medium  was  mixed  for  30  minutes  and  then  allowed  to  settle  for  5  minutes.  The 
supernatant  was  then  separated  for  placement  in  eight  123  mL-reactors.  Each  of  the 
eight  reactors  was  fed  with  6  mL  of  the  supernatant  and  topped  with  buffered  dilution 
water  to  bring  up  the  final  volume  to  60  mL.  The  control  reactor  was  topped  up  to  62 
mL.  All  the  reactors  were  then  capped  with  potassium  hydroxide  pellets  in  the  holders 
attached  to  the  caps. 

The  capped  reactors  were  then  placed  in  the  respirometer  water  bath,  and  oxygen 
was  supplied  to  each  reactor  for  a  period  of  4  hours.  During  this  period  the 
temperature  of  the  water  bath  was  maintained  constant  at  25‘’C,  and  the  contents  of 
the  reactors  were  kept  well-mixed  by  magnetic  stirrer  bars.  At  the  end  of  the  4-hr 
period,  each  reactor,  except  the  control,  was  administered  with  different  volumes  of  the 
toxicant  being  assayed.  Each  toxicant  dose  was  prepared  by  dissolving  the  toxicant  in 
2  mL  of  acetone.  The  dosed  reactors  and  the  control  reactor  were  then  placed  in  the 
respirometer,  with  continued  supply  of  oxygen.  The  data  acquisition  system  was  then 
re-initiated  to  monitor  the  oxygen  uptake  of  each  reactor  for  the  next  6  hour  period. 

For  each  toxicant,  a  trial  run  was  first  made  to  establish  the  range  of 
concentrations  to  bring  about  50%  inhibition.  Once  this  range  was  determined  by  trial 
and  error,  a  final  run  was  made  with  appropriate  volumes  of  the  toxicant  to  cover 
inhibitions  in  the  range  of  about  30  to  80%.  The  oxygen  uptake  as  a  function  of  time  as 
recorded  by  the  data  acquisition  system  was  then  used  to  determine  the  IC5Q  values 
as  described  in  the  main  text. 


Appendix  III 

Detailed  Statitics  of  OSAR  Models 


Regression  Summary  Residual  Statistics 

log  Y-RM  vs.  1XV  _  log  Y-RM  vs.  1XV 


Count 

14^ 

#>-0 

7 

Num.  Missing 

0 

#<0 

7 

R 

SS[e(i)-e(l-1)l 

1.665 

R  Squared 

.910 

Ourbin-Watson 

2.288 

Adjusted  R  Squared 

.903 

Serial  Autocorrelation 

-.236 

RMS  Residual 

.246 

ANOVA  Table 
log  Y-RM  vs.  1XV 


DF  Sum  of  Squares  Mean  Square  F-Value  P-Value 


Regression 

1 

7.398 

7.398 

122.014 

<.0001 

Residual 

12 

.728 

.061 

Total 

13 

8.126 

Regression  Coefficients 
log  Y-RM  vs.  1XV 


Coefficient 

Std.  Error 

Std.  Coeff. 

t-Value 

P-Value 

Intercept 

3.690 

206 

3.690 

17.920 

<.0001 

1XV 

-.896 

.081 

-.954 

-11.046 

<•0001 

Confidence  Intervals 
log  Y-RM  vs.  1XV 

Coelltciem 

Intercept  ' 

1XV 


95%  Lower  95%  Upper 


3.690 

3241 

4.138 

-.896 

-1.072 

-.719 

Regression  Summary 
log  Y-RM  vs.  1XV 

Split  By:  Type 

Cell:  Alkane 

Residual  Statistics 
log  Y-RM  vs.  1XV 

Split  By:  Type 

Cell:  Alkane 

Count 

5 

#>-0 

3 

Num.  Missing 

0 

#<0 

2 

R 

.999 

SS[e(i)-e(i-l)l 

.003 

R  Squared 

.999 

Durbin-Watson 

2578 

Adjusted  R  Squared 

QQQ 

Serial  Autocorrelation 

-.403 

RMS  Residual 

.018 

ANOVA  Table 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Alkane 


Regression 

Residual 

Total 


DF  Sum  of  Squares  Mean  Square  F-Value  P-Value 


Regression  Coefficients 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Alkane 


Coefficient  Std.  Error  Std.  Coeff.  t-Value  P-Value 


Intercept 

1551 

.044 

1.851 

41.741 

<-0001 

1XV  _ 

-.765 

.014 

-.999 

-53.133 

<.0001 

Confidence  Intervals 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Alkane 


Regression  Summary 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Amine 
Count 

Num.  Missing 
R 

R  Squared 
Adjusted  R  Squared 
RMS  Residual 


6 

0 

.957 

.915 

.894 

.101 


Residual  Statistics 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Amine 
#>-0 
#<0 

SS[e(i)-e(i-l)l 
Ourbin-Watson 
Serial  Autocorrelation 


3 

3 

.104 

2.549 

-.326 


ANOVA  Table 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Amine 


DF 

Sum  of  Squares 

Mean  Square 

F-Value 

P-Value 

Regression 

1 

.443 

.443 

43.327 

.0028 

Residual 

4 

.041 

.010 

Total 

5 

.484 

Regression  Coefficients 
log  Y-RM  vs.  1 XV 
Spin  By:  Type 
Cell:  Amine 


Coefficient 

Std.  Error 

Std.  Coeff. 

t-Value 

P-Value 

Intercept 

1.045 

.135 

1.045 

7.738 

.0015 

1XV 

-.470 

.071 

-.957 

-6.582 

.0028 

Confidence  Intervals 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Amine 


Intercept 

1XV 


Coefficient  95%  Lower  95%  Upper 


1.045 

.670 

1.420 

-.470 

-.669 

-.272 

Regression  Summary 
log  Y-RM  vs.  1XV 
Spin  By:  Type 
Cell:  Aromatic 


Residual  Statistics 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Aromatic 


Count 

9 

#>-0 

5 

Num.  Missing 

0 

#<0 

4 

R 

.852 

SSIe(i)-e(i-1)l 

1.125 

R  Squared 

.726 

Durbin-Watson 

1.665 

Adjusted  R  Squared 
RMS  Residual 

.687 

.311 

Serial  Autocorrelation 

.053 

ANOVA  Table 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Aromatic 


DF 

Sum  of  Squares 

Mean  Square 

F-Vakje 

P-Value 

Regression 

1 

1.792 

1.792 

18566 

.0035 

Residual 

7 

.676 

.097 

Total 

8 

2.468 

Regression  Coefficients 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Aromatic 


Coefficient 

Std.  Error 

Std.  Coeff. 

t-Value 

P-Value 

Intercept 

3.258 

.737 

3.258 

4.420 

.0031 

1XV 

-1.133 

.263 

-.852 

-4509 

.0035 

Confidence  Intervals 
log  Y-RM  vs.  1XV 
Split  By:  Type 
Cell:  Aromatic 


Interce... 

1XV 


Coefficient  95%  Lo...  95%  Upper 


3258 

1515 

5.001 

-1.133 

-1.755 

-.511 

Regression  Summary 
log  Y-RM  2  vs.  OXV 

Count 

Num.  Missing 
R 

R  Squared 
Adjusted  R  Squared 
RMS  Residual 


12 

39 

.942 

.887 

.876 

.141 


Residual  Statistics 
log  Y'RM  2  vs.  OXV 
#>.0 
#<0 

SSIe(i)-e(i-1)] 
Duitoin-Watson 
Serial  Autocorrelation 


6 

6 

.394 

1.990 

-.066 


ANOVA  Table 
log  Y-RM  2  vs.  OXV 


DF 

Sum  of  Squares 

Mean  Square 

F-Value 

P-Value 

Regression 

1 

1558 

1558 

78.803 

<0001 

Residual 

10 

.198 

.020 

Total 

11 

1.756 

Regression  Coefficients 
log  Y-RM  2  vs.  OXV 


Coefficient 

Std.  Error 

Std.  Coeff. 

t-Value 

P-Value 

Intercept 

2.670 

239 

2.670 

11.178 

<.0001 

OXV 

-.448 

.050 

-.942 

-8B77 

<•0001 

Confidence  Intervals 
log  Y-RM  2  vs.  OXV 

^  Coefficient 
Intercept  " 

OXV 


95%  Lower  95%  Upper 


2.670 

2.137 

3.202 

-.448 

-.560 

-335 

Appendix  IV  a 

Results  ot  10-Component  Mixtures 
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Appendix  IV  b 

Results  of  8-Component  Mixtures 
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Appendix  V 

Results  of  8-Component  Mixtures 
Comparison  of  Experimental  Results  vs.  Predictions 


Appendix  V-  Results  of  8-Component  Mixture  Studies 
Comparison  of  Experimental  vs.  Predicted  Concentrations  Causing  50%  Inhibition 
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